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ABSTRACT 

We have determined the cluster X-ray temperature function from two flux- and redshift-limited samples 
of clusters. The first sample is comprised of 25 clusters with average redshift 0.05. The local temperature 
function derived from it supercedes the one we previously published (Henry & Arnaud, 1971). Fourteen 
clusters with average redshift 0.38 comprises the second sample. We perform maximum likelihood fits of 
cluster evolution models to these data in order to constrain cosmological parameters. For an open model 
with zero cosmological constant we find that the density parameter is Q = 0.49 ± 0.12, the rms mass 
density fluctuation averaged over 8 h 1 Mpc spheres is a 8 = 0.72 ± 0.10, and the effective index of the 
mass density fluctuation spectrum on cluster scales is n = -(1.72 ± 0.34) where all errors are symmetrized 
at 68% confidence. The corresponding results for the case where a cosmological constant produces a flat 
universe are: 0.44 ± 0.12, 0.77 + 0.15, and -(1.68 ± 0.38). These results agree with those determined from 
a variety of different independent methods, including supernovae, galaxy-galaxy correlations, fluctuations 
in the microwave background, gravitational lens statistics, and cluster peculiar velocities. 

Subject headings: cosmology: observations-galaxies: clusters: general-large-scale structure of universe- 
X-rays: galaxies 

1. INTRODUCTION 

The properties of clusters of galaxies are rather sensitive to several parameters of cosmological interest. 
One of the first to be constrained was a 8 , the amplitude of mass density fluctuations on a scale of 8 h" 1 
Mpc. Here, h is the current value of the Hubble parameter H in units of 100 km s 1 Mpc" 1 . Using the 
number density of nearby clusters as a function of the velocity dispersion of their member galaxies, 
Evrard (1989) found a value between 0.67 and 1, Cole & Kaiser (1989) found 0.5 ± 0.1 and Peebles, Daly 
& Juszkiewicz (1989) found ~ 0.67. Henry and Arnaud (1991, HA91) were the first to use cluster X-ray 
temperatures to measure a 8 , finding 0.59 ± 0.02 from a local sample. 

For reasons of simplicity, not to mention fashion, these early analyses were restricted to Q = 1, where Q 
is the ratio of the mass density of the universe at the present epoch to the critical density required to halt 
the expansion. More recent work has relaxed that assumption, showing that a local sample can only 
constrain a combination of Q ar »d <J 8 . A large number of papers have performed these kinds of analyses 
using cluster abundances as a function of velocity dispersion, temperature, or mass (deduced from either 
the previous two observables or from gravitational lensing), or the logN - logS number counts (White, 
Efstathiou, & Frenk 1993; Eke, Cole, & Frenk, 1996; Viana & Liddle 1996; Kitayama & Suto, 1997; 
Mathiesen & Evrard, 1998 among others). The result of all this work is that a 8 is in the range of 0.5 to 
0.6, with errors ranging from 5% to 20% at 68% confidence, if Q = 1- As Q becomes less than one a 8 
increases, in order that the number of peaks above a given threshold remains at the observed value. It is 
worth noting that the value a 8 is reasonably robust, having remained close to the above range since the 
very first determinations. 



Additional information is required to break the degeneracy between a 8 and Q - O ne source of that new 
information is the cosmic microwave background fluctuations. These fluctuations are currently measured 
on a vastly different size scale from clusters and some method has to be used to bridge the gap with its 
associated error. Another source of additional information is the evolution of cluster number density, 
which has the advantage that the new data pertain to the same size scale as the old. It was realized early 
on that cluster evolution was extremely sensitive to Q (White and Rees, 1978; Perrenod, 1980; Evrard, 
1989; Peebles et al. 1989; Frenk et al. 1990; Lilje, 1992; Oukbir & Blanchard 1992). The few high 
redshift clusters known, as well as instrumental difficulties, limited the early application of this technique 
using X-ray observations (Perrenod & Henry, 1981; Henry et al, 1982). Cluster evolution is not as 
sensitive to the cosmological constant, parameterized by A,o, the present ratio of the cosmological constant 
A to 3 times the square of H (Eke et al, 1996; Viana & Liddle, 1996). 

Recently, a number of papers have appeared that constrain Q through a comparison of the properties of 
high and low redshift clusters (Carlberg et al, 1997; Bahcall, Fan & Cen 1997; Fan, Bahcall & Cen, 1997; 
Henry, 1997, H97; Sadat, Blanchard & Oukbir, 1998; Eke et al, 1998; Viana & Liddle, 1999; Borgani et 
al, 1999; Reichart et al, 1999; Donahue & Voit, 1999; Kay & Bower 1999). There is a wider divergence 
of results compared to the situation with <j 8 . With the exception of Sadat et al (1998) and Reichart et al 
(1999), all this work is consistent with a low density universe (0.2 < £2 < 0.5) with Q = 1 excluded at 
varying levels of significance ranging from less than one sigma to greater than 10" 6 ! Some of this 
discrepancy is surely due to the different redshift ranges probed; higher redshift objects have more power 
to discriminate. Almost just as surely, however, are uncorrected systematic biases, either observational or 
theoretical. 

The shape of the mass fluctuation spectrum may be measured as well. Early work showed that the 
standard cold dark matter model was inconsistent with the cluster distribution, either as a function of 
temperature (HA91; Oukbir & Blanchard, 1997), X-ray luminosity and its evolution (Henry et al, 1992) 
or of mass derived from temperatures, and/or velocity dispersions (Bahcall & Cen 1992). Rather a steeper 
spectrum (more negative spectral index) was preferred. Constraining the shape of the fluctuation spectrum 
has not received as much attention as the other cosmological parameters, with most workers assuming a 
cold dark matter spectrum with shape parameter measured from galaxy surveys. However, cluster data 
provide independent constraints that agree with those from galaxies. 

We have previously given a brief report (H97) of the constraints on the above cosmological parameters 
provided by the cluster X-ray temperature function and its evolution. Using the same data, we have also 
investigated the effects of systematic uncertainities on these constraints (Eke et al, 1998). In this paper we 
improve on H97 as well as provide all the details of its analysis. Among the improvements are new X-ray 
temperatures from ASCA for the low redshift sample, and an enlarged high redshift sample for which the 
temperatures of all clusters have been derived uniformly from the same analysis using the latest ASCA 
calibrations. We also compare our results with those from independent constraints. 

2. OBSERVATIONAL RESULTS 

Our primary observation goal is to construct the cluster temperature functions at low and high redshifts. 
Since rich clusters have a wide range of temperatures, from ~2 keV to -15 keV, we can only measure any 
evolution by comparing the temperature distributions. There are no standard temperature clusters. 

The observed temperature function is a quite steep function of temperature (approximately [kT] 5 ) and it is 
difficult to show it graphically with small sized samples. After experimenting with several methods, we 
have decided to present our data in the form of the integral temperature function, n(>kT). This form has 
the advantage that it is nonparametric and unbinned. Its disadvantage is that the error on any point is 
difficult to establish since all points but the hottest are correlated. However, because of the very steep 
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dependence on temperature, even a factor of two difference in temperature is enough to insure virtually 
uncorrelated data Thus we show representative error bars calculated from the Poisson noise on the 
number of clusters used to determine a given point on the distribution. The previous remarks refer to the 
graphical representation of our results. In this paper we will always use a maximum likelihood fitting 
technique to the individual unbinned measurements that properly calculates the errors on the fit 
parameters. 

The integral temperature function is estimated from the data using 

n{>m=Y,^— a) 

i " sea J 

where the sum is over all clusters with kTj > kT and V sea>i is the search or maximum volume in which the 
i th cluster could have been detected given its selection criteria. 

2.1 LOW REDSHIFT SAMPLE 

The low redshift sample remains the 25 clusters from HA91. We give the properties of these objects in 
Table 1. This sample came originally from the early all-sky surveys with non-imaging detectors. 
Confusion has been reduced by using pointed observations, still with non-imaging detectors, but with 
greatly reduced fields of view compared to the survey instruments. All but the Coma and Virgo clusters 
are unresolved with these instruments, so accurate values of the total flux are available. 

The selection criteria are that the total flux in the 2 - 10 keV band is greater than 3 x 10 11 erg cm 2 s" 1 , the 
redshift is less than 0.17, and the absolute galactic latitude is greater than 20° (plus a small zone of 
exclusion around the Magellanic Clouds) yielding a solid angle surveyed of 8.23 steradians. The redshift 
restriction was designed to insure a local sample. Its precise value was set to exclude clusters with 
temperatures greater than 10 keV through the luminosity - temperature relation (see Section 3). The 
maximum redshift at which a cluster could have been detected is given to an adequate approximation by 
the Euclidean limit 
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where F lim is 3 x 10 11 erg cm 2 s Use of the Euclidean expression for the search volume calculated from 
z max introduces an error of less than 30%. Still we use the exact formula for an open model given in the 
appendix for the low redshift sample. This calculation is an improvement on our previous work. The 
difference between the search volumes for open and flat models is about 5% for this sample, which we 
neglect. We give the search volumes in Table 1 and the resulting integral temperature functions in the 
table and Figures 1 and 2. The median search volume for the low redshift sample is 2.9 x 10 7 h" 3 Mpc 3 . 

We need to mention the case of A2163. This cluster satisfies all the selection criteria except its redshift is 
0.201, thus it is excluded from our low redshift sample. Certainly, it would be difficult to argue that a 

cluster at such a redshift is local. At a temperature of 14.6^ keV it would have been the hottest in our 

sample. Its z max is 0.213 implying a search volume of 4.96 x 10 8 h 3 Mpc 3 (for Q = 0.5), so if it were 
included in the sample, it would have added 2.02 x 10" 9 h 3 Mpc" 3 to all points in the integral temperature 
function. This contribution is small, except for the hottest cluster in Table 1 . 
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An important improvement over H97 is the recent availability of temperatures from ASCA measurements 
of nearly two-thirds of the sample (Markevitch et al 1998). We want to compare temperature functions at 
low and high redshifts, so we use the ASCA single temperature fits uncorrected for cooling flows from 
Markevitch et al (1998) because such corrections are not yet possible for the high redshift clusters. Until 
now the temperatures for the low and high redshift samples came from different instruments. These new 
measurements help to eliminate a possible source of systematic error. Fortunately for previous work, the 
ASCA and non-ASCA temperature scales are the same to a few percent, as is expected. The error- 
weighted mean ratio of ASCA to mean non-ASCA temperatures in Table 1 for clusters that can be 
compared is 1.031 ± 0.022 (68% confidence error on the mean). 

Edge et al (1990), HA91, and Markevitch (1998) have previously reported the low redshift temperature 
function. Their differential temperature functions were fit to a power law, n(kT) = A[kT] a . The best 
fitting A (in units of h 3 Mpc 3 keV" 1 ) and a were: &.&™xl0~* and 4.93+0.37, 1.8^jjcl0" 3 and 4.7±0.5, 

and 3.5!2 2*10~ 4 and 4.2+0.7 respectively. The power law indices agree within their errors. The 

normalizations are somewhat difficult to compare because of the steep temperature dependence. At a 
temperature of 6 keV, the differential temperature functions are (in units of 10 7 h 3 Mpc 3 keV" 1 ) 

1.28+0.08, 3.96^ 3 g2> and 1.90+0.54 respectively after allowing for the correlated errors in A and a (a 

larger a requires a larger A to maintain the total number of clusters). 

This comparison shows that the HA91 normalization is systematically high, although only marginally so 
with respect to the Markevitch (1998) determination. Eke et al (1996) pointed out to us that this higher 
normalization was in error. Upon investigation, we found that the HA91 differential temperature function 
contained two nearly compensating errors. The first was a typographical error in a numerical constant in 
our code. The second was subtler (almost anything is subtler than a typo). A determination of the 
differential temperature function requires binning the data and performing an average over the objects in 
the bin. This procedure will introduce some arbitrariness that the use of the integral function avoids. 
There are two obvious averages that come to mind. The one used in HA91 was l/<V seai AT;> where <. . .> 
denotes average. The estimator commonly used is <l/(V seaji AT ; )>. The former average is always less than 
the latter, unless all objects in the average are identical. Further, the contribution of a given cluster to the 
temperature function depends on that of other clusters with similar temperatures using the first average. In 
the second average, the contribution of any clusters to the temperature function depends only on its 
unique properties, which is clearly to be preferred. 

We have, therefore, made a new fit of the low redshift temperature function to a power law using the 

maximum likelihood procedure described in section 4. The best fitting parameters are A = xl0~ 3 

h 3 Mpc 3 keV 1 and a = 4.9+0.7. This function is practically indistinguishable from that of Edge et al 
(1990) and Markevitch (1998). 

2.2 HIGH REDSHIFT SAMPLE 

The high redshift sample comes from the Einstein Extended Medium Sensitivity Survey (EMSS, Gioia et 
al 1990; Henry et al 1992, H92 hereafter; Gioia and Luppino 1994). The EMSS was conducted with an 
imaging detector; hence its fundamental variable is surface brightness. Sources in the EMSS were found 
serendipitously in fields targeted at unrelated objects. Each source was required to produce a sufficient 
number of counts in a 2.4' x 2.4' detect cell to be significantly above the expected background in that cell 
for the exposure of the field. The selection criteria to be included in the sample used here are that the 
detect cell flux in the 0.3 - 3.5 keV band be greater than 2.5 x 10 13 erg cm" 2 s" 1 , 5(1950) > -40°, and the 
redshift satisfy 0.3 < z < 0.6. There are sixteen clusters in Gioia and Luppino (1994) satisfying our 
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selection criteria. However two of them, MS 1333.3 and MS 1610.4, are X-ray point sources (Molikawa 
et al, 1998; our unpublished analysis). We removed these two objects from our sample leaving a total of 
fourteen. All objects discussed here are extended X-ray sources in deep ROSAT HRI observations, that is 
they are all bona fide clusters. Finally we have revised the redshift of MS 1241. 5 based on new spectral 
observations taken at the Hawaii 2.2m telescope (Gioia, 1998 private communication). After these 
revision the high redshift sample has an average redshift of 0.38. We also give our measurements for 
MS2053.7 whose flux is just below the cutoff and thus not in the sample used to constrain cosmological 
parameters. 

We have determined temperatures, fluxes, and luminosities for the high redshift clusters from ASCA 
observations. These data are given in Table 2. All clusters were analyzed identically using the latest 
calibrations. Specifically, we have used version 2 of the telescope effective area and point spread 
functions. The temperatures were hardly affected compared to H97 and Eke et al (1998). However, these 
files eliminate the ad hoc 18% correction to the Gas Imaging Spectrometer (GIS) flux employed there, 
although comparison of the two results shows that it was an excellent fudge factor. We have also taken 
care to use the 128 channel GIS detector 3 stuck bits Redistribution Matrix File (RMF) for MS 0451.6, 
the only cluster in our sample so affected. Spectra were extracted from 6. 125' and 2.5' radii regions for the 
GIS and Solid-state Imaging Spectrometer (SIS) respectively. Background was acquired from adjacent 
regions of the same exposure. We grouped the spectra until each bin contained at least 20 counts, and then 
fit the two GIS and two SIS spectra simultaneously using the program XSPEC 10. The adjustable 
parameters were the GIS and SIS normalizations, temperature, abundance, and Hydrogen column density 
for the SIS (because of low energy calibration uncertainties). Fixed parameters were the redshift and GIS 
Hydrogen column density. The total flux and luminosity were determined from the GIS normalizations, 
since most of the observations were conducted in 2CCD mode and thus some of the photons were lost to 
gaps between the chips and off the edge of the detectors. We assumed that these objects are point sources 
in order to determine a total flux, a good approximation for distant clusters observed with the GIS on 
ASCA. The results are in Table 2. Our measurements of the temperature are in excellent agreement with 
previous values reported by others. The fluxes determined by ASCA in the 0.3 - 3.5 keV band are -15% 
higher here than in H97 and Eke et al (1998) because they mistakenly did not correct for interstellar 
absorption. The 2-10 keV luminosities are hardly effected by this omission. 

The standard EMSS procedure for determining the fraction, f, of a cluster's total flux in the EMSS detect 
cell is to integrate a beta model cluster surface brightness distribution, with [3 = 2/3, over the detect cell 
and divide that result by the flux from the same beta model integrated to infinity, see equation (1) of H92: 



where D is the angular half-size of the detect cell, 1.2', and o (z) is the angular size of the cluster core 
radius, a ( /D A (Sl >z) with D A the angular diameter distance, see the Appendix, and a the linear core radius. 
We use a = 0.125h 1 Mpc, which is appropriate for low (Mohr, Mathiesen & Evrard, 1999) and high 
(Vikhlinin et al, 1998) redshifts. The error weighted mean ratio of the ASCA to EMSS flux for all sources 
in Table 2 is 1.17 + 0.05 (68% confidence error on the mean). The median value is 1.19. These values 
apply to either total or detect cell fluxes, since the same f from equation (3) converts either the EMSS 
detect cell to total or the ASCA total to detect cell. The factor ~2 discrepancy for MS0015.9, the largest 
deviation, was first pointed out by Nichol et al (1997) and may result from the object lying at the extreme 
edge of the field of view in the EMSS. Excluding this source gives a weighted average flux ratio of 1.12 + 
0.05 (68% confidence error on the mean) and a median of 1.15. We show a comparison of the fluxes of 
these clusters from the original EMSS and our new ASCA measurements in Figure 3. 




(3) 
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The cause of this small discrepancy is not yet clear. We can think of a number of causes. First, it may 
simply be due to the absolute calibration uncertainties of the Einstein and ASCA observatories which are 
of the same order. Second, it could result from unrelated point sources in the GIS extraction beam that has 
a solid angle 20 times larger than the EMSS detect cell. To assess this possibility, we have examined the 
HRI images of all sources in Table 2. There are five objects (MS2137.3, MS1241.5, MS1358.4, 
MS0353.6, and MS 1224.7) that do not have obvious point sources in the GIS extraction beam. All others 
have one or more such sources. The error weighted mean ratio of the ASCA to EMSS flux for these five 
is 0.99 ± 0.08 (68% confidence error on the mean). The median value is 1.09. Therefore, we regard not 
being able to exclude unrelated sources in the much larger GIS beam as the likely cause of the 
discrepancy. 

However, Jones et al (1998) and Ebeling et al (1999) note a somewhat larger discrepancy in the same 
sense between the EMSS and ROSAT PSPC and HRI determined fluxes respectively and so we 
investigated further. One effect not included in equation (3) is the IPC point spread function. Including 
the effects of the point spread function on f, f PSF , is substantially more complicated. In the end, f/f P sF is 
nearly a constant at all redshifts. For z > 0.3, it is 1.373 to better than 0.3%; for all redshifts for which the 
EMSS was used for evolutionary studies, it deviates from 1.373 by no more than 4.5%. This correction 
would explain the discrepancies noted above, particularly the larger ones reported by Jones et al (1998) 
and Ebeling etal (1999). 

Once we have opened the Pandora box of trying to determine a value of f that is accurate to -10% other 
effects must also be included. No cluster has X-ray emission at infinite radius, as we have assumed so far. 
In fact, the emission will only extend to the virial radius, given according to the theory outlined in Section 
4 by, r v (z v ,kT) = (3M/[47iA(n ,z v ) £l Q p c , (l + z v ) 3 ]) 1/3 = 1.90 h" 1 Mpc (3 TM 1/2 (kT/10keV) 1/2 (1 + z v ) 3/2 (for 
Q = 1) = 2.08 h" 1 Mpc (kT/10keV) 1/2 (1 + z v )" 3/2 (for |3tm given in Section 4). Most of these terms are 
defined after equation (12). For now we note that z v is the virialization redshift that we take to be the 
redshift of the cluster. In addition, Vikhlinin, Forman, & Jones (1999) have shown that the surface 
brightness in the outer regions of clusters falls off faster than that given by our standard (3 = 2/3. This 
value was determined with Einstein Imaging Proportional Counter (IPC) measurements that were not as 
sensitive as the ROSAT Position Sensitive Proportional Counter (PSPC) that is able to trace the cluster 
emission out to larger radii. From the data in Vikhlinin et al (1999) we find (3(kT) = (0.621 + 0.022) + 
(0.017 ± 0.004) kT. Including the PSF, the new value of |3 and retaining the fixed core radius given above, 
but only integrating out to the virial radius yields a new f, f P sFkT, that is -6% larger than f when averaged 
over the ranges 0.3 < z < 0.6 and 4 < kT < 10 appropriate to our situation. The range of f PS FkT over the 
same range of z and kT is 0.98 - 1.13. While a correction of this size would probably resolve the EMSS - 
ASCA flux discrepancy, it is likely too small for that noted by Jones et al (1998) and Ebeling et al (1999). 
And still we have not yet included the different core radii of different clusters. Furthermore, for our 
purposes, we only need the one-half power of the total flux to determine the volume accessible to a given 
temperature cluster (see Section 3.2), which gives only a few percent change from the standard EMSS f. 
Given the remaining uncertainties and the smallness of the changes after including everything described 
above, we have retained the original EMSS fraction given by equation (3). 

Given the wide range of exposures on the targets in the EMSS, there is a correspondingly wide range of 
solid angles surveyed to a particular flux limit. There is not a single flux limit for the entire survey region 
as there is for the low redshift sample. The integral sky coverage, £l surv , as a function of limiting flux, is in 
Table 3 of H92. We repeat it here in our Table 3 and also give the differential sky coverage, d£2 surv> that is 
the additional solid angle surveyed upon going to the next higher flux limit. The maximum redshift, z max , 
to which a cluster at redshift z and detect cell flux F det can be detected for the i th flux limit, is given by 
equation (3)ofH92: 
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(4) 



where D L (£2 ,z) is the luminosity distance (see the Appendix). Note that in equations (3) and (4) above Q 
= 1 is to be used because the value of a was derived using that assumption. The volume searched for each 
cluster is then obtained by summing over all flux limits: 



sea 



dV(Q ,<rrun(z M ,z max ,)) dV(Q < Zl ) 



da 



(5) 



where z h and z u are the lower and upper redshift boundaries respectively (0.3 and 0.6 for our case) and 
dV(Q. ,<z)/dQ. is the volume per unit solid angle within a redshift z given in the Appendix. 

We give these volumes in Table 2 using the detect cell flux derived by multiplying the total flux from our 
ASCA observations by the f appropriate for each cluster. By using the ASCA data to determine the detect 
cell flux instead of the original EMSS data, we have decided to accept a -15% systematic uncertainity 
instead of -25% statistical uncertainity. The volumes were calculated assuming Q = 0.5. The median 
search volume for clusters in the high redshift sample is 6.4 x 10 7 and 7.7 x 10 7 h 3 Mpc 3 for a zero 
cosmological constant universe (termed open) and for a spatially flat universe in which Q. + X Q = 1 
(termed flat) respectively. These volumes are comparable to that for the low redshift sample. We also give 
the integral temperature functions in Table 2 and in Figures 1 and 2 as calculated from equation (1). 

3. MAXIMUM LIKELIHOOD METHOD AND TEMPERATURE SELECTION FUNCTION 

3.1 Maximum Likelihood 

We use a maximum likelihood fit to the unbinned data in order to determine various model parameters. 
The method is described in Marshall et al (1983). The likelihood function is given by their equation (2), 
which for our situation is: 



S =-2£ln 



n{Q. ,Zi,kTi) 



rfV(Q ,z,-) 
dzd£l 



kT max 



+ 2 \dkT J" n(Q ,z,kT)a(z,kT) 



kTmin 



d 2 v{a ,z) 

dzdQ. 



dz (6) 



Here N is the number of clusters observed, n(£2 ,z,kT) i s the temperature function, £2(z,kT) is the solid 
angle in which a cluster with temperature kT at redshift z could have been detected (the selection 
function) and d 2 V/dzd£2 is the differential volume, which is given in the Appendix. The d 2 V/dzd£2 in the 
first term of equation (6) is not in Marshall et al (1983) because it is a function of one of the parameters 
(£l ) we are trying to constrain. The usual procedure is to assume a value for this parameter. This term 
was not in the likelihood function of H97. Its inclusion has only a small effect; the best fit values with and 
without it agree within their 68% confidence errors. However, we retain it here. 

The best estimates of the model parameters and derived quantities are obtained by minimizing S. 
Confidence regions for them are obtained by noting that S is distributed as % 2 with the number of degrees 
of freedom equal to the number of interesting parameters (Lampton, Margon and Bowyer, 1976; Avni 
1976). We have either one, two or three interesting parameters depending on what we are trying to 
constrain. Sometimes we are only interested in £l , at other times we want to know the full temperature 
function or to derive constraints on Q ar >d CJ 8 (three parameters). 
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The maximum likelihood fit only provides the best fit, it does not provide an assessment of whether that 
fit is a good fit to the data. Marshall et al (1983) require that the two independent cumulants 



C(< z,) = 



rkT max 
'kT min 



n(Q , z, kT)Cl(z, kT)(d 2 V (Q , z) I dzdQ.)dkT 



(7) 
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n(Q , z, kT)£l(z, kT)(d 2 V(£l , z) I dzdQ.)dkT 



'kT min 



and 



C(> kT i \z i ) = 



rkl max / r\ 

\ m n(Ll ,z 



rkT max //-v 




.,kT)£l( Zi ,kT)dkT 
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both be uniform according to the two-tailed Kolmogorov-Smirnov test. They recommend rejecting the 
model if the probability of either being uniform is less than 0.05. 



The clusters in both the low and high redshift samples were selected from flux limited surveys. They were 
not selected from a temperature-limited survey; such a survey may not even be possible. Fortunately, 
clusters obey a temperature - luminosity relation and the luminosity may be derived from the redshift and 
flux. Thus we recast the selection as a function flux into a selection as a function of redshift and 
luminosity and, via the temperature - luminosity relation, as a function of redshift and temperature. 

There have been many discussions of the cluster temperature - luminosity relation. Among the larger 
and/or more recent compilations are David et al. (1993, 104 clusters), White, Jones, and Forman (1997, 
86 clusters), Mushotzky and Scharf (1997, 38 clusters), Markevitch (1998, 35 clusters), Allen and Fabian 
(1998, 30 clusters), Arnaud and Evrard (1999, 24 clusters) and Henry and Sornig (1999, 29 clusters). The 
relation from the data in Table 2 is shown in Figure 4. The conclusions from all this work are that clusters 
have a definite temperature - luminosity relation that exhibits little or no evolution at any redshift but the 
relation has a large scatter related to the size of the cooling flow at the cluster center. Not all of these 
compilations are applicable to our work here because some have included the effects of cooling flows in 
their analysis, either by excluding the central region, or by choosing samples that contain mostly cooling 
flow or non-cooling flow clusters. It is not yet generally possible to discern the existence of cooling flows 
at high redshift, much less correct for them. 

We give in Table 4 a selection of recent results that are most relevant to this paper. The two low redshift 
samples analyze essentially all available data and use mostly the same clusters. The Henry and Sornig 
(1999) result comes from all publicly available ASCA observations of clusters as of July 1, 1998 with 
redshifts between 0.3 and 0.6 and contains the objects in this paper as well as others. None of the fits 
given differentiate between cooling flow and non-cooling flow clusters. It is apparent that there is little 
evolution between z = 0.39 and 0.06 and therefore we use the David et al (1993) fit to determine our 
selection functions for both samples. 

There are two additional complications. First, in order approximately to compensate for the large scatter 
in the temperature - luminosity relation we have used the -1 a values of both the normalization and slope 
of the David et al (1993) fit. Second, we convert the fit in Table 4 from the bolometric luminosity to the 
band in which the clusters were selected. Specifically, for the low redshift sample we use kT = 
3.55[L44(2,10)] 0286 and for the high redshift sample we use kT = 3.05[L44(0.3,3.5)] a355 , after allowing for 
both effects. The selection functions are given in Figure 5, using f in equation (3) for the high redshift 
case. A comparison with the data from Tables 1 and 2 indicates they provide a plausible description of the 
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selection. The most discrepant cluster is C0336, one of only two non-Abell clusters in the low redshift 
sample. This object is an extended X-ray source not contaminated by point sources (Irwin & Sarazin, 
1995). 

4. THEORETICAL TEMPERATURE FUNCTION 

The primordial fluctuation spectrum is the squared modulus of the spatial Fourier transform of the 
fluctuations in the underlying matter density field about their average. We use a power law model for this 
quantity: 



P(k)= ' 



(9) 



In this expression the shape parameter is n and V u is a large volume in which the statistical properties of 
the fluctuations are the same as in any other such volume. The normalization is set by ko. Nearly always 
the normalization is parameterized with the rms mass fluctuations on an 8h 1 Mpc scale, a 8 . This power 
spectrum gives the rms value of the density fluctuations at zero redshift when smoothed with a top hat 
containing a mass M: 



aJa o ,0,M) = 0.675 



r(3 + n)sin(n7t/2) 
2 n n(2 + n)(l-n)(3-n)_ 



p.5k h~ l Mpc) Q„ hM l5 \ (10) 



(Peebles, 1980) where M 15 is the mass in units of 10 15 solar masses. Then, a 8 =a p (£2 o ,0,0.594 x 10 15 h l Q. D ) 
which is independent of Q . 

The spatial distribution of the fluctuations that grow by gravity into the structure we see is determined by 
P(k). The rate of their growth is determined mainly by £2 . What objects can we use to trace this 
evolution? Probably the best are X-ray emitting clusters of galaxies. The physics and astrophysics of X- 
ray emission from clusters of galaxies are relatively straightforward. Simple gravitational processes 
dominate cluster formation and evolution and imply that clusters are still forming today. A cluster itself is 
non-linear, having collapsed, but it formed out of a background mass density field that is still evolving 
linearly, an enormous simplification. These same formation processes also heat gas trapped by the cluster 
potential, which then produces optically thin thermal radiation. The gas is in simple coronal equilibrium; 
that is the ion populations are mostly in the ground state and electron impact dominates the ionization 
balance and excitation processes. 

Given the simplicity of the relevant processes, there is hope that the evolution of cluster X-ray emission 
can be reliably calculated and those calculations may be verified by direct observations of nearby objects. 
This situation is in stark contrast to that for supernovae, quasars and galaxies, the only other objects that 
are observable at cosmological distances. The evolution of these objects is too complicated to be reliably 
calculated and/or they formed too long ago for it to be well observed. 

The evolution of the cluster mass function is given by the Press-Schechter (1974, PS) formula. Although 
relatively simple, PS has been repeatedly shown to provide a very accurate description of N-body 
calculations (e.g. Lacey and Cole, 1994; Mo and White, 1996), at least over the range of masses of 
interest to us (Gross et al, 1998). The number of objects per unit mass interval per comoving volume is, 
according to PS: 



9 



n(Cl ,z,M) = - 




2Q p c , 8 c (Q ,z) D(Q o ,Q) JlnG p (Q o ,0,M) 
71 M o p (Q o ,0,M) D(Cl ,z) dM 



exp 



-6 2 (O o ,z)D 2 (C2 o ,0) 



(11) 



2c 2 p (Q o ,0,M)D 2 (n o ,z) 



where p c o is the present critical density, 8 c (£2 ,z) is the mass density fluctuation required for a spherical 
perturbation to collapse at redshift z, and D(£\,,z) is the growth factor of perturbations while they are 
evolving linearly. Recent calculations have now derived these quantities for most cases of interest. We 
give them in the Appendix for the three model universes we consider here. 

However, the cluster mass is not an easily measurable quantity, although it may be determined from X- 
ray observations in some circumstances. More readily observable are the X-ray luminosity or 
temperature. The X-ray luminosity depends on the square of the hot gas density, which may have its own 
evolution independent of that of the cluster. Perhaps the best cluster observable readily accessible to 
current instrumentation is the integrated temperature. As one might suspect, and numerical hydrodynamic 
calculations confirm the suspicion, the average cluster temperature is given by the depth of the potential 
well of the dark matter of the cluster. If the gas were hotter (cooler) than that given by virial equilibrium, 
then it would have blown away (collapsed) long ago because the time scale for changes in the 
configuration of the gas, the sound crossing time, is much shorter than the Hubble time. 

The temperature function is obtained from the mass function of equation (11) via the chain rule; n(kT) = 
n(M) dM/dkT. The mass-temperature relation is given by considering a modified singular isothermal 
sphere: 



Here A(£2 Q ,z v ) is the ratio of the cluster's average density to that of the background mass density at the 
virialization redshift (see the Appendix) and [3 TM is the modification factor accounting for departures from 
virial equilibrium. Bryan & Norman (1998, their Figure 5) show that equation (12) is a good description 
of the mass - temperature relation in numerical hydrodynamic simulations over at least the redshift range 



We determine |3 TM from numerical hydrodynamic simulations of clusters, the results from which we give 
in Table 5. The average value is 1.21 + 0.05. In an intriguing paper on this issue, Shapiro, Iliev and Raga 
(1998) have derived a value of 1.16 by considering a truncated isothermal sphere that has minimum 
energy. Physically, the pressure to provide the truncation comes from the outward going shock front in a 
spherical infall solution. This new work provides the final piece of a completely analytic model of cluster 
evolution. In the future it may be desirable to determine [3 TM empirically, but at present such 
determinations have a quite substantial error as the last three rows in Table 5 show. For the moment, the 
empirical determinations do not disagree with those from the hydrodynamic simulations or the analytic 
result. 

Combining equations (10), (11) and (12) gives the evolving cluster temperature function. 



kT = 



1.9Skev Q A(Q ,Q 



(hM l5 ) 2n (l + z v ). 



(12) 



0-1. 
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n(£l , z, kT) = 1 .29xl(T 5 (h^Mpc)' 3 (keVy 1 Q S C (Q , z) f ( ^° ,Q) (1 + z v ) (3 "" )/4 (9.5h- y MpcQT^k () ) (3+ " )/2 



rp™^r )/4 



7.98 



£>(a o ,z) 

x(n-3)/12 ( t 1/2 

18ti 2 I |2"n(l-n)(3-n)(2 + n)(3 + n) 2 I 



[Q A(Q O J 



exp 



•i.i5 2 (a ,z) 



£> 2 ayo) 

D 2 (£l ,z) 



(1+Zv) 



-(3+n)/2 



r(3 + n)sin(«7l/2) J 
^ R 7 ^ A ( 3+ »"V 187C 2 



7.98 



x(3+n)/6 



^ A(a ,zj J 



(9.5/* Mpc^ fc ) . — — 

r(3 + n)sin(n7i/2) 



(13) 



We will also need the integral function, which is 



n(Cl ,z,>kT) = 2.58*1 (r^/T'Mpc)-' (1.1) 



,(3-n)/(2«+6) 



7.98 



P 



TM 



8 e (O .z) 



£>(n ,z) 



6/(n+3) 



(9.5h~ l Mpc k ) 3 



2"n(l-n)(3-n)(2 + n) 
r(3 + n)sin(mt/2) 



3/(n+3) 



H — 3 

2(n + 3) 



,1.1 



5 c (a ,z) 



£>qyo) 



(l+Zv) 



-(b+3)/2 



18ti 2 



n(3+«)/6 



(9.5h- l MpcQ.~ U3 k ) 



3+n f2 n n(l-«)(3-n)(2 + n)]f $ m kT T +V ' 2 



r(3 + /z)sin(mc/2) 



7.98 



(14) 



We adopt the late collapse approximation and set the redshift at which the cluster is observed equal to the 
virialization redshift, z = z v . It is possible to relax this approximation, but the effects of doing so are not 
large. Figure 8a of Kitayama & Suto (1996), Figure 5 of Kitayama & Suto (1997) and Figure 3 of Viana 
and Liddle (1999) compare various results with and without the approximation. Note that the parameter s 
= in both papers by Kitayama & Suto corresponds to our equation (12). In general, there is almost no 
change at all for either the flat model (for redshifts under consideration here) or the open model for large 
Q.q. For example, the inferred value of Q f° r the flat model changes by ~ 0.01 near Q = 0.3 and z = 0.4, 
judging from Figure 3 of Viana and Liddle (1999). Only for Q ~ 0.3 open models at z ~ 0.4 is the effect 
on Q ~0.1, but this size is still comparable to the statistical uncertainity of the present samples. 

5. MEASURING COSMOLOGICAL PARAMETERS 



We are now ready to fit the data presented in Section 2 to the theory described in Section 4 using the 
maximum liklihood method of Section 3. The best fitting integral temperature functions are shown in 
Figures 1 and 2. The two cumulants given in equations (7) and (8) are uniform with a probability of 0.70 
and 0.14 respectively for the open model and 0.78 and 0.13 respectively for the flat model. Thus each fit 
is acceptable. 
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The most general description of the results requires the three parameters of the fit. We show these results 
in Figures 6 and 7. It is straightforward to read off the value of n, which is -(1.72^3') and 

-(1.68!q4i) at the 68% confidence for the open and flat models respectively. These values agree well 
with that measured in HA91. The shape parameter that we adopt, n, may be converted to the shape 
parameter of the popular cold dark matter fluctuation spectrum, T, using equation (5) of Pen (1998) or 
equations (2) and (3) of Viana and Liddle (1999). Both give nearly the same results, which are that our fit 
implies T « 0.05 ± 0.22. 

The presentation in Figures 6 and 7 is somewhat difficult to appreciate, so we also give the constraints for 
fewer parameters. First we concentrate just on £l in Figure 8. Constraints are relatively tight when 
considering this single parameter. We find that Q = 0.49 ^° f at the 68% confidence level and at the 

95% confidence level for the open model. For the flat model the corresponding values are Q = 0.44 ^j'jj 

at the 68% confidence level and at the 95% confidence level. These results agree with those in H97, 

except the errors are ~ 40% smaller here. A closed universe is ruled out at much greater than 99% 
confidence. 

We recast the constraints in Figures 6 and 7 into a more conventional format in Figures 9 and 10. Three 
parameters are still required, but the constraints on n and ko are collapsed into a 8 . Some of the degeneracy 
between these two parameters still remains in the form of a banana shaped contour particularly for the flat 

model. Still, we are able to make a relatively precise measurement of a 8 , finding values of 0.72^,J )8 and 

0.77 ^'iq at 68% confidence level for the open and flat cases respectively. These results also agree with 
those in H97. 

6. COMPARISON WITH OTHER MEASUREMENTS 

Determining the values of cosmological parameters is one of the fundamental goals of astrophysics. 
Recently, the number of results has increased dramatically and it is impossible to provide a 
comprehensive review in one section of a paper reporting still more. We have given in the Introduction a 
short summary of what was obtained using clusters of galaxies utilizing similar methods to ours. Here we 
want to give a few of the most recent results using methods that are very different. 

Our analysis in terms of the fluctuation spectrum, P(k), permits a direct measurement of this quantity at 
spatial frequencies approximately 0.2 to 0.5 h Mpc The positions of galaxies, either as a two 
dimensional angular distribution or a three dimensional spatial distribution, give a direct measurement of 
the shape (assuming the bias between galaxies and mass is independent of scale, eg Mann, Peacock, & 
Heavens, 1998), but not the normalization, of P(k) over a much wider range. Gaztanaga & Baugh (1998) 
have measured the local slope of P(k) using the angular correlation function of the APM Galaxy Survey. 
Galaxies are in the non-linear regime at the spatial frequencies probed by clusters, but after correcting for 
this effect the APM measurements imply that n = -(1.6 ± 0.2), which is in excellent agreement with our 
result. 

The baryon fraction contained in clusters, combined with £2 b provided by Big Bang nucleosynthesis, has 
been used to measure the same three cosmological parameters that we have. Evrard (1997) finds Q = 
(0.30 + 0.07) h 2/3 , which is Q = 0.42 ± 0.10 for h = 0.6. Shimasaku (1997) determined the cluster gas 
mass function from which he derived (o 8 ,n) = (0.64_ 009 ,— 1.6^ 04 ) for h = 0.5 and 

(o 8 , n) = (0.73!°; fi -1 -5!o 4 ) for h = 0.8. All of these results agree well with ours. 
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Deviations from a smooth Hubble flow, termed redshift distortions or peculiar velocities, caused by the 
gravity of mass fluctuations constrain the quantity £2 CT 8 /o 8 (tracer) , where Q%{tracer) is defined 

analogously to a 8 using the objects tracing the mass field and velocities in it. The degeneracy between the 
parameters measured by this method is similar to that from cluster abundances at a single redshift. The 
peculiar velocity constraint has been applied mostly to samples of galaxies. Three recent analysis of the 
IRAS 1.2 Jy sample, where a 8 (IRAS) = 0.69 ± 0.04 (Fisher et al, 1994), are described by Sigad et al 
(1998), da Costa et al (1998), and Willeck et al (1998). They find O.° 6 O & = 0.61 ± 0.09, 0.41 ± 0.07 and 
0.34 ± 0.05 respectively. The discrepancy among these three analyses of nearly the same data is 
somewhat unexpected. A similar measurement using clusters of galaxies gives ^o' 6 O g = 0.44^ 'j' (90% 

confidence, Watkins, 1997; Borgani et al, 1997). We compare the da Costa et al (1998) result with ours in 
Figures 11 and 12. 

The statistics of the number of gravitationally lensed quasars can constrain the cosmological constant. 
Some recent work and its 95% confidence limit (unless otherwise noted) on Q , assuming the flat model, 
is by Falco, Kochanek & Munoz (1998; > 0.26), Cooray (1999; > 0.21), Cooray, Quashnock & Miller 

(1999; > 0.21), Helbig et al (1999; > 0.36), Chiba & Yoshii (1999; 03™ at 68% confidence) and Cheng 
& Krauss (1999; 0.25 - 0.55 systematic effects dominate). All these results agree with what we find. 

Fluctuations in the cosmic microwave background (CMB) determined with COBE normalize P(k). This 
normalization may be extrapolated from the scales probed by COBE to cluster scales if the shape of P(k) 
is known. Bunn & White (1997), among many others, use a CDM spectrum to perform the extrapolation. 
Their results are shown in Figures 11 and 12. The CMB data nicely complement that from clusters and 
both combine to yield a small overlap region for either the open or flat models that results in a somewhat 
lower value of Q than our best fit. 

A classical test used to determine cosmological parameters is the Hubble diagram. Collins and Mann 

(1998) present a K^band Hubble diagram for brightest cluster galaxies in X-ray luminous clusters that has 
very small scatter. They find good agreement with us, with H = 0.28 ± 0.24 for the open model and 
0.55^5 for the flat model (both at 68% confidence). 

Of course one of the most productive uses of the Hubble diagram recently has come from the one 
constructed from Type la supernovae. Efstathiou et al (1999) combine this diagram with measurements of 
fluctuations in the CMB at COBE and smaller scales. They exclude the open model at high confidence 
but the flat model is quite consistent with their data. If we limit their results to such a model, in order to 
compare with our own, they find Q = 0.29 + 0.08 (at 68% confidence). Tegmark (1999) has emphasized 
that including additional physics, such as gravity waves, will substantially increase the size of the error 
bars. Nevertheless, these results are consistent with ours. 

Although the above has not been comprehensive, we hope it has been representative. Most recent 
analyses indicate that £2 lies approximately in the range 0.25 - 0.55. As useful as the comparison of 
results among different methods is, a better approach is to combine several data sets in order to determine 
a joint best fit. This approach has at least three advantages: smaller errors, breaking of degeneracies 
among the parameters, and possible revelation of systematic errors. Webster et al (1998), Lineweaver 

(1999) , and Bridle et al (1999) have taken initial steps in this direction. The latter work performs a 
simultaneous fit to the CMB fluctuations, the IRAS 1.2 Jy redshift survey and the data presented here. 
There is good agreement among all three data sets and with our results. 

7. FUTURE WORK 
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All but one EMSS cluster with z > 0.3 has now been observed with ASCA. The last one is unobservable 
with ASCA due to its proximity to NGC4151, which is very bright. Thus it may be possible to repeat the 
analysis here with ~ 20 high redshift clusters, a number comparable to the present low redshift sample. In 
addition to improving the statistics because of increased sample size, some of the additional clusters will 
be at higher redshifts than any used here and that provides increased leverage on the measurement. 

Another statistical error is that on the temperature measurement. Since the temperature function is a very 
steeply falling function, simple statistical errors on temperature bias it high. Further, the errors on 
temperature are larger for the high redshift sample compared to the low, which may lead to some bias. We 
are probably close to the Time Allocation Committee limit, so it is not feasible to improve greatly on this 
situation with ASCA. However, the upcoming XMM mission will provide a large increase in effective 
area that should substantially improve the statistical quality of the temperatures. 

Useful as an improvement in statistics will be, we believe that the present measurement is starting to be 
limited by systematic uncertainties. We have already discussed in Section 4 the effects of assuming all 
clusters collapse at the epoch they are observed. Another systematic effect comes from the cluster 
temperature - luminosity relation, which is integral to the determination of the selection function. This 
relation has a very large scatter, but the scatter is reduced considerably if the effects of the central cooling 
flow are removed (Allen & Fabian, 1998; Markevitch, 1998). It is only possible to make this exclusion 
with ASCA for clusters at low redshift, unless the observation has a sufficiently large number of detected 
photons that multiple component spectral models can be fit. As mentioned before, the ASCA exposures 
are already quite long so this approach is not possible for us. However, both the Chandra and upcoming 
XMM missions will permit a simple spatial exclusion of the cooling flow. Robinson, Gawiser, & Silk 
(1999) and Eke et al (1998) suggest another method of limiting the systematic effects caused by the 
temperature - luminosity relation. The hottest clusters in a redshift limited sample are detectable in the 
entire volume defined by the redshift cuts because they are so luminous. Thus the temperature - 
luminosity relation is mostly irrelvant. Viana & Liddle (1999) also discussed the utility of working only 
with the hottest clusters. Unfortunately, adopting this approach with the present samples yields very large 
statistical errors because there are so few clusters hot enough. 

We believe the best way to eliminate the effects of systematic errors is to repeat the analysis performed 
here with completely new samples. Several such samples, both at low and high redshift, have now been 
compiled from ROSAT. A relatively modest investment of XMM time, comparable to what has already 
been used here with ASCA, would permit a totally independent measurement of the same cosmological 
constants to better statistical precision. 

Notwithstanding the desirable additional work described above, we want to end on a positive note. Most 
independent methods imply that the universe will expand forever. We have an empirical determination of 
the mass fluctuation spectrum. Clusters abundances and CMB fluctuations provide the normalization at 
the small and large scales and galaxy redshift surveys provide the shape in between. The combination of 
these two tells the story of how the magnificent large scale structures we see about us came to be. 

I am pleased to acknowledge many useful discussions with Harald Ebeling and Vincent Eke, which 
greatly improved my understanding of the maximum likelihood technique, the EMSS selection function 
and the theory of cluster evolution. E. Miyata is the Co-PI of the MS1621.5 and MS2053.7 observations. I 
also have benefited from conversations with K. Arnaud, M. Arnaud, N. Bahcall, A. Blanchard, S. Cole, S. 
De Grandi, M. Donahue, A. Evrard, C. Frenk, I. Gioia, M. Hattori, J. Irwin, A. Liddle, C. Mullis, J. 
Navarro, M. Norman, J. Ostriker, J. Primack, R. Sadat, P. Shapiro, Y. Suto, M. Turner, P. Viana, 
A.Vikhlinin, M. Voit, and S. White. Over the years this work has been supported by NASA through 
grants NAG5-2523 and NAG5-4828. 
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APPENDIX 



Here we collect the various quantities, defined in the body of the paper, that are required for our 
theoretical analysis of the temperature function and maximum likelihood fits to our data. As a reminder, 
8 c (£2 ,z) is the critical overdensity for collapse at redshift z; A(Q Q ,z) is the ratio of the average cluster 
mass density to that of the background density for a cluster which collapses at redshift z; D(£2 Q ,z) is the 
growth factor up to redshift z; D L (£2 D ,z) is the luminosity distance to redshift z and the angular diameter 
distance is D L (£2 ,z)/(l+z) 2 ; V is the volume and dCl is the solid angle surveyed. Quantities with subscript 
zero are evaluated at the present epoch. We distinguish three cosmologies, closed, open and flat. 

Closed Model £l = 1 



5 C (U) = ^^~ 1.686 
20 



(Peebles, 1980, Eq. 19.48a evaluated at 2k) 



A(l,z) = 187T 2 =177.7 (Peebles, 1980, Eq. 19.50 multiplied by 32. A 

factor of 8 results from assuming the cluster collapses to a radius half its value at maximum expansion 
and another factor of 4 from the decreasing background density during the time required to collapse that 
is twice the time to reach maximum radius.) 



D(l,z) = 



1 



l + z 



(Peebles, 1980, Eqs. 11.5,11.7 and 92.2) 



D L {\,z) = 2— (1 + z-VTTI) 



dV(l,<z) = 8 
da 3 



( \ 
c 

V J 



1- 



VI 



+ z 



d 2 V(l,z) 
dzdQ. 



= 4 



' c Y 2+z-2VT+z 



V J 



(l + z) 



5/2 



Open Model Cl < 1, A, = 



„2^ 
l + z 

2 



fl(z) = 



2 + x 



15 



5 C (£l ,z) = -(l + (2ti) 2/3 {jx 2 + 2x - In [(1 + x) + Vx 2 +2x \ 
2 1 



2/3 > 
J 



( 



1 + 



6 _3V2 + xln[(l + x) + Vx 2 + 2x 



3/2 



A(a ,z) = 



(27i) 2 x 3 

{\/x 2 +2x - In (l + x) + Vx 2 +2xjf 
1993 Eq A15; Oukbir and Blanchard, 1997) 



v 



(Maoz, 1990 Eq 12; Lacey and Cole, 



2x n 



. 6 3^2 + x\r\l + x) + V x + 2x 

1 + — 175 

X X 



(Peebles, 1980 Eq. 11.16; our x is twice 



his; the different normalizing factor yields the same form as £l — > 1 as the closed model.) 
2 c 



D L (a ,z) = - Y ^-[a z + (2-a )(i-^i+a z)] 



a 2 h 



dV(£l <z) 



c 

v J 



1 



2(1 -Q ) 



3/2 



In 



Vi+q z + a /i-q 



Vi+a z-Vi-a J [i-Vi-^o , 



In 



+ 



2-a (2-a +a z)Vr+iv 



Q 2 (i-C2 ) a 2 (i-a )(i+z)^ 



+ l u ( 2 [(2-a )(Vi+a z-i)-a zj 

" (1+Z) 



This equation is numerically unstable at £l a = 1 and for low £2 D and/or z. Consequently, it is more reliable 
to integrate the following equation over z. 



d 2 V(£l ,z) 
dzda 



= 4 



f A 3 
c 



[q z + (2-a )(i- A /i+^ z)] 2 
cit(i+z)^i+n i)Z 



Flat Model ft + X = 1 

\l/3 



X = 



_(a -'-i) 1 



G(z) = 



1+z 

1 

1 + x 3 
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5 C (Q , z) = 3(12 ^ )2/3 [l - 0.0123 log(l + x' )] 
Eke et al, 1996) 



A(Q , z) = 187t z (1 + 0.4093x z,n ' 2 ) 
their w is our x 3 ; Eke et al, 1996) 



D(a ,z) = 

Eq. 14) 

D(a ,z) = 



X 



Jl^7Cdy(l + x 3 y 6,5 y 3 < 2 

JO 



(Nakamura and Suto, 1997, Eq. C-28; 



(Nakamura and Suto, 1997, Eq. C-19 



(Peebles, 1984, 
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cos -1 
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(Einsenstein, 1998, Eqs. 8 - 10) 



where K = sin 75° = 0.5V2 + V3 . Here F(cp,K) and E(cp,K) are incomplete Legendre elliptic integrals of 
the first and second kind. Sometimes the following recursion relation must be used when cp > n/2: F(cp,K) 
= 2 F(Jl/2,K) - F(7l - cp,K) (Abramowitz and Stegun, 1964, equations 17.3.2 and 17.4.3) and similarly for E. 



dy 



D L (a Q ,z) = —(l+z) , 

H J p (j (l+y) 3 +(l-a ) 

Eq. 12; Caroll, Press, & Turner, 1992, Eqs. 23 and 25) 



(Yoshii, 1995, 
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Eq. 14) 

The integral in the last three equations is 
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JJ 



(Yoshii, 1995, 



(see Mathiesen and Evrard, 1998 and references in Eisenstein, 1998). The numerical evaluation of the 
integrals is probably somewhat simpler to code than the Legendre closed form. Feige (1992) also gives an 
extensive discussion of the use of elliptical integrals in cosmological formulas. 
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We graph the critical overdensity for collapse in Figure 13. Note that this parameter is never more than a 
few percent different from its closed value (Lilje, 1992). We show in Figure 14 the growth factor to 
redshift z relative to the present epoch. It is only this combination that enters into the Press-Schechter 
mass function. Here, different values of the density parameter result in a large difference in growth. This 
difference is amplified by the exponential term of the Press-Schechter mass function. Finally, we show 
the cosmology dependent factors of the mass-temperature relation in Figure 15. 
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FIGURE CAPTIONS 



Figure 1. Integral cluster temperature function with best fitting open model superposed. The open (filled) 
circles are the low (high) redshift cluster observations, while the dotted (solid) lines are the low (high) 
redshift best-fitting models, respectively. 

Figure 2. Integral cluster temperature function with best fitting flat model superposed. The open (filled) 
circles are the low (high) redshift cluster observations, while the dotted (solid) lines are the low (high) 
redshift best-fitting models, respectively. 

Figure 3. Ratio of the ASCA to EMSS total fluxes as a function of the total flux measured with ASCA. 
The open circles are for clusters with no obvious point sources in their HRI images within 6.125' of the 
cluster center. 

Figure 4. The data with ± la errors are the luminosity-temperature relation for clusters in the redshift 
range 0.3 - 0.6 reported here. The X-ray luminosity is in the 2-10 keV band and is calculated using h = 
0.5 and £2 = 1 to be consistent with previous low-redshift work. The two solid lines are the ± la fits 
reported by David et al. (1993) for their low-redshift sample without any shift in normalization. 

Figure 5. The selection functions for the low- and high-redshift cluster samples in the redshift- 
temperature plane with the data from Tables 1 and 2 superposed. The dotted lines on the high redshift 
selection function correspond to the EMSS integral sky coverage in Table 3, with the smallest solid angle 
at the coolest temperature. The solid angle for the low redshift selection function is a constant 8.23 
steradians. 

Figure 6. Sixty-eight percent confidence contours for the three interesting parameters n, k) and Q , 
corresponding to an increase in liklihood of 3.5 from the minimum, for the open model. Also shown are 
lines of constant a 8 . 

Figure 7. Sixty-eight percent confidence contours for the three interesting parameters n, k) and Q > 
corresponding to an increase in liklihood of 3.5 from the minimum, for the flat model. Also shown are 
lines of constant a 8 . 

Figure 8. Delta liklihood contours as a function of the single interesting parameter Q - The dotted lines 
show various confidence levels. 

Figure 9. Sixty-eight percent confidence contour for a 8 and Q for the open model. The constraints on n 
and k) have been collapsed into a 8 , so three free parameters must still be considered. The contour thus 
corresponds to an increase in liklihood of 3.5 from the minimum. 

Figure 10. Sixty-eight percent confidence contour for a 8 and £2 for the flat model. The constraints on n 
and k) have been collapsed into a 8 , so three free parameters must still be considered. The contour thus 
corresponds to an increase in liklihood of 3.5 from the minimum. 

Figure 1 1 . The comparison of the constraints provided by the data presented in this paper with those from 
galaxy redshift distortions and the COBE normalization of a cold dark matter fluctuation spectrum for the 
open model. 

Figure 12. The comparison of the constraints provided by the data presented in this paper with those from 
galaxy redshift distortions and the COBE normalization of a cold dark matter fluctuation spectrum for the 
flat model. 
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Figure 13. The critical overdensity for collapse as a function of the present value of the density parameter 
and redshift for the open and flat models. 

Figure 14. The linear theory growth factor from a redshift z to the present as a function of the present 
value of the density parameter for the open and flat models. 

Figure 15. The factors dependent on cosmology in the mass-temperature relation as a function of the 
present value of the density parameter and redshift for the open and flat models. 
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Name 

A2142 

A2029 

A401 

A1656 

A754 

A2256 

A399 

A3571 

A478 

A3667 

A3266 

A1651 

A85 

A119 



0.0894 
0.0766 
0.0736 
0.0232 
0.0543 
0.0581 
0.0723 
0.0397 
0.0882 
0.0556 
0.0590 
0.0846 
0.0556 
0.0442 



F(2,10) a 

7.5 
7.5 
5.9 
32.0 
8.5 
5.2 
3.4 
11.5 
6.6 
6.7 
5.9 
3.7 
6.4 
3.0 



EXOSAT 

kT(keV) 

11.0±1.9 



8.0+0.3 
8.7±1.7 



7.6+1.1 
6.8+1.1 
5.3+1.6 



TABLE 1 
LOW REDSHIFT SAMPLE 

Einstein GINGA ASCA 

kT(keV) kT(keV) kT(keV) 

9.4+1.0 8.68+0.20 8.8+0.6 

7.8+1.2 8.7+0.3 

7.8+1.0 8.3+0.5 

8.3+0.9 8.21+0.09 

9.1+1.1 7.57+0.31 9.0+0.5 

7.3+0.7 7.51+0.19 7.5+0.4 

5.8+1.3 7.4+0.7 

10.9+2.9 6.9+0.3 

7.8+1.0 6.70+0.19 7.1+0.4 



5.1+0.9 



6.5+1.6 
6.2+0.6 

6.2+0.5 
5.9+1.0 



7.0+0.6 
7.7+0.8 
6.3+0.5 
6.1+0.2 
5.8+0.6 



AVE 
kT(keV) 
8.74+0.19 
8.65+0.29 
8.20+0.45 
8.19+0.09 
8.04+0.25 
7.50+0.17 
7.04+0.62 
6.99+0.29 
6.80+0.17 
6.76+0.53 
6.74+0.48 
6.30+0.50 
6.11+0.19 
5.65+0.45 



Open (flo = 0.5) 

V sea b n(>kT) c 

1.63E8 0.00E+0 

1.06E8 6.13E-9 

6.77E7 1.56E-8 

2.8 1E7 3.03E-8 

4.80E7 6.59E+8 

2.89E7 8.68E-8 

2.94E7 1.21E-7 

3.02E7 1.55E-7 

1.32E8 1.88E-7 

3.66E7 1.96E-7 

3.61E7 2.23E-7 

5.20E7 2.51E-7 

3.42E7 2.70E-7 

5.90E6 3.00E-7 



REF 

15,6,5,5,11 

19,5,11 

3,5,11 

3,6,5,10 

3,6,5,2,11 

3,5,2,11 

3,5,11 

13,6,9,11 

19,6,5,18,11 

3,14,5,11 

3,5,11 

19,11 

3,5,11 

3,6,5,11 
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A3558 


0.0479 


4.2 


3.8+2.0 




5.70+0.20 


5.5+0.3 


5.63+0.17 


1.22E7 


4.69E-7 


3,5,4,11 


A1795 


0.0630 


5.3 


5.1±0.5 


5.8+0.5 


5.34+0.11 


6.0+0.3 


5.42+0.10 


3.73E7 


5.51E-7 


3,6,5,5,11 


A2199 


0.0299 


7.1 


4.7+0.4 


4.5+0.3 


4.54+0.10 




4.54+0.09 


6.63E6 


5.78E-7 


3,6,5,1 


A2147 


0.0353 


3.3 


4.4+1.7 


4.4+0.4 






4.40+0.39 


3.50E6 


7.29E-7 


15,6,5 


A3562 


0.0502 


3.5 


3.8+0.9 


4.2+1.2 






3.94+0.72 


1.07E7 


1.01E-6 


13,6,9 


A496 


0.0328 


5.7 


4.7+0.9 


3.9+0.2 


3.91+0.06 




3.91+0.06 


6.30E6 


1.11E-6 


3,6,5,18 


A3526 


0.0103 


11.2 


3.6+0.4 


3.9+0.2 


3.54+0.13 




3.64+0.11 


5.63E5 


1.27E-6 


3,6,5,18 


A1367 


0.0213 


3.5 


3.5+0.5 


3.7+0.3 


3.50+0.18 




3.56+0.14 


8.64E5 


3.04E-6 


3,6,5,1 


A1060 


0.0124 


4.4 


3.3+0.4 


3.9+0.4 


3.27+0.05 


3.1+0.2 


3.27+0.05 


2.43E5 


4.20E-6 


3,6,5,8,17 


C0336 


0.0349 


4.7 


3.1+0.3 


3.0+0.1 


3.45+0.10 




3.22+0.07 


5.70E6 


8.32E-6 


12,6,9,16 


Virgo 


0.0039 


30.0 


2.4+0.3 


2.6+0.1 


2.20+0.02 




2.21+0.02 


1.35E5 


8.49E-6 


7,6,9,2 



a 10" 11 erg cm 2 s 1 in the 2 to 10 keV band 
b h" 3 Mpc 3 
c h 3 Mpc" 3 

All fluxes are from Edge, Stewart, & Fabian (1992), all temperature errors are at the 90% confidence level. 

References: 1. Arnaud (1994); 2 Arnaud & Evrard (1998); 3 den Hartog & Katgert (1996); 4 Day et al (1991); 5 David et al (1993); 6 Edge 
& Stewart (1991); 7 Federspiel, Tammann, & Sandage (1998); 8 Hatsukade (1989); 9 HA91; 10 Hughes et al (1993); 11 Markevitch et al 
(1998); 12 McNamara, Bregman, & O'Connell (1990); 13 Melnick & Moles (1987); 14 Piro & Fusco-Femiano (1988); 15 Struble & Rood 
(1987); 16 Tsuru (1993); 17 Tamura et al (1996); 18 Yamashita (1992); 19 Zabludoff, Huchra, & Geller (1990) 



25 



TABLE 2 
HIGH REDSHIFT SAMPLE 

F det (0.3,3.5) L 44 (2,10) a 

(10 13 erg cm 2 s" 1 ) ^ erg s" 1 kT Open (Q = 0.5) Flat (Q = 0.5) 

MS Z EMSS ASCA ASCA (keV) V sea b n(>kT) c V sea b n(>kT) c 

0451.6 0.5392 9.5+1.7 12.4±0.2 30.6+0.6 in V 09 L41E8 0.00E+0 1.80E8 0.00E+0 

0015.9 0.5466 7.1+0.8 14.1+0.2 34.1+0.6 8.9+0.6 1.42E8 7.09E-9 1.81E8 5.57E-9 

1008.1 0.3062 5.9+0.8 9.2+0.3 8.7+0.3 8 9 +12 6.86E7 1.41E-8 8.29E7 1.11E-8 

°- z '-i.i 

1358.4 0.3290 12.2+2.1 11.0+0.3 10.2+0.3 6.9+0.5 9.13E7 2.87E-8 1.14E8 2.32E-8 

1621.5 0.4274 3.4+0.6 5.2+0.2 7.5+0.3 fifi*' 9 7.11E7 3.96E-8 8.96E7 3.20E-8 

u - u -0.8 

0353.6 0.320 6.3+1.3 7.8+0.3 6.8+0.3 6 5 +L0 5.83E7 5.37E-8 7.09E7 4.31E-8 

1426.4 0.320 4.4+0.6 5.3+0.2 4.6+0.2 fi 4 +1 -° 2.69E7 7.09E-8 3.24E7 5.72E-8 

1241.5 0.549 4.2+1.1 4.7+0.3 10.7+0.6 6 l +1 - 4 1.02E8 1.08E-7 1.30E8 8.81E-8 
1147.3 0.303 3.0+0.7 5.3+0.2 4.1+0.2 f. n +1 ° 2.13E7 1.18E-7 2.55E7 9.58E-8 
0811.5 0.312 2.6+0.5 3.5+0.2 2.6+0.2 4 a+ 10 6.98E6 1.65E-7 8.34E6 1.35E-7 

2137.3 0.313 19.3+2.6 20.9+0.4 15.3+0.3 4.9+0.3 1.27E8 3.08E-7 1.60E8 2.55E-7 

0302.7 0.4246 3.8+0.4 4.2+0.3 4.8+0.3 4 4 +0 - 8 4.89E7 3.16E-7 6.10E7 2.61E-7 

1224.7 0.3255 5.3+1.1 4.2+0.2 3.0+0.2 A l +07 1.64E7 3.36E-7 1.99E7 2.78E-7 

^- 1 -0.5 

1512.4 0.3726 4.5+1.0 4.5+0.3 3.8+0.3 3.4+0.4 3.47E7 3.97E-7 4.26E7 3.28E-7 
2053.7 0.583 2.48+0.58 2.4+0.2 6.6+0.6 8 1 +37 

°- 1 -2.2 

a D. = 1 b h 3 Mpc 3 c h 3 Mpc 3 

All errors are at the 68% confidence level. Redshifts with four figures are from Carlberg et al (1996), 
those with three figures are from Gioia & Luppino (1994) except MS 124 1.5 is from Gioia (private 
communication). 
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TABLE 3 

HIGH REDSHIFT SAMPLE SKY COVERAGE 



F det (0.3,3-5) Differential Solid Angle (dQ surv ) Integral Solid Angle (£l surv ) 

(lCT 13 erg crn 2 s" 1 ) (ster) (ster) 

2.50 0.0551 0.0551 

2.78 0.0097 0.0648 

3.34 0.0182 0.0830 

4.01 0.0234 0.1064 

4.81 0.0281 0.1345 

5.77 0.0267 0.1613 

6.93 0.0232 0.1845 

8.31 0.0177 0.2022 

10.0 0.0104 0.2126 

11.9 0.0062 0.2188 

14.3 0.00295 0.2218 

17.3 0.00128 0.2230 

20.7 5.5 x 10 4 0.2236 

24.8 1.5 x 10~ 4 0.2237 
29.7 3 x 10 5 0.2238 
35.7 3 x 10" 5 0.2238 
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TABLE 4 

TEMPERATURE - BOLOMETRIC LUMINOSITY RELATION 2 

Reference <z> A a kT @ L 44 (bol) = 7 h Clusters 

David et al (1993) 0.066 2.94+0.14 0.297±0.004 5.24+0.20 104 

White, Jones & Forman (1997) 0.057 2.76+0.08 0.33+0.01 5.25+0.05 86 

Henry & Sornig (1999) 0.385 3.23+0.45 0.253+0.046 5.28^ 29 

a kT = A [L 44 (bol)] a for H = 50 km s" 1 Mpc 1 , ft = 1, and kT in keV 

b errors allow for correlated errors on A and a, ie -1 a on A is associated with +1 a on a and vice-versa. 
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TABLE 5 



DETERMINATION OF THE MASS - TEMPERATURE PARAMETER (3 



Method 










h 




Objects 


Ref 


Hydrodynamic Simulation 


1 26+0 03 


1.0 


0.10 





0.50 


0.59 


22 


1 


MvHroHvnamip Si mill a ti on 


1 17 

X , X / 


1 

X . V7 


10 


o 


50 


59 63 


42 


2 


1-TvHroHvnamir Simulation 


0.88 


0.2 


0.10 


0, 0.8 


0.50 


1.00 


16 


2 


HvdroHvnflmic Simulation 


i i c+n nc 


1.0 


0.05 


o 


0.50 


0.50 


20 


3 


1-TvHrodvnjimip Simulation 

X iy Ul UUJ llcll IX1C* LJXXXXU.Xdl.Xl_/Xl 


1 07+0 04 


3 


04 


7 


70 


1 05 


30 


4 


Hvdrodvnamic Simulation 


1.29 


1.0 


0.06 


o 


0.50 


1.05 


25 


5 


rxjLilVJLljllallil^ ollllUiaLlVJll 


1 36 


1 


10 


o 


fiS 


67 


95 




Hydrodynamic Simulation 


1.36 


1.0 


0.075 





0.50 


0.70 


25 


5 


Hydrodynamic Simulation 


1.31 


0.4 


0.06 





0.65 


0.75 


25 


5 


Truncated Isothermal Sphere 


1.16 


1.0 













6 


Gravitational Lens 


1.02±0.11 


1.0 







0.50 




8 


7 


Velocity Dispersion 


1.21±0.35 








1.0 
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8 


Spatially Resolved Temperature 










1.0 




11 


8 



References: 1. Evrard (1990); 2. Evrard, Metzler, & Navarro (1996); 3. Pen (1998); 4. Eke, Navarro, & 
Frenk (1998); 5. Bryan & Norman (1998); 6. Shapiro, Iliev, & Raga (1998); 7. Hjorth, Oukbir, and van 
Kampen (1998); 8. Horner, Mushotzky, & Scharf (1999) 
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